Universal scheme to generate metal-insulator transition in disordered systems 
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We propose a scheme to generate metal-insulator transition in random binary layer (RBL) model, 
which is constructed by randomly assigning two types of layers. Based on a tight-binding Hamilto- 
nian, the localization length is calculated for a variety of RBLs with different cross section geometries 
by using the transfer-matrix method. Both analytical and numerical results show that a band of 
extended states could appear in the RBLs and the systems behave as metals by properly tuning the 
model parameters, due to the existence of a completely ordered subband, leading to a metal-insulator 
transition in parameter space. Furthermore, the extended states are irrespective of the diagonal and 
off-diagonal disorder strengths. Our results can be generalized to two- and three-dimensional disor- 
dered systems with arbitrary layer structures, and may be realized in Bose-Einstein condensates. 

PACS numbers: 71.30.+h, 71.23.An, 72.20.Ee, 03.75.-b 



One of the most important issues in condensed-matter 
physics is Anderson localization 1], which predicts that 
the electronic wavefunctions may become localized in im- 
perfect crystals and leads to a disorder-induced metal- 
insulator transition (MIT), owing to the quantum inter- 
ference between multiple scatterings of an electron with 
random impurities and defects 2-4]. Another seminal 
work along this direction is the scaling theory of local- 
ization which indicates that all electronic states are 
exponentially localized in low-dimensional noninteracting 
systems even for infinitesimal disorder and become local- 
ized in three-dimensional (3D) systems with sufficiently 
large disorder strength. In fact, Anderson localization 
is a universal wave phenomenon and has been observed 
experimentally in a wide variety of systems, including 
light 0-0]; microwaves [ioj . acoustic waves 11, 12 1, and 



matter waves [131-41 7 1 



The simplest one-dimensional (ID) example exhibiting 
Anderson localization is the random binary alloy, where 
the two constituents A and B are randomly distributed in 
the lattice [see Fig.QJa)] and all eigenmodes are localized 
within a small region [see the solid line in Fig.QJc)] . How- 
ever, Anderson localization breaks down in the random 
binary alloy when correlations are introduced in the dis- 
order distribution. A discrete number of extended states 
have been reported in the random-dimer model 18- 2fjj| 



and its generalized versions [2l|, |22j , where one or both 
sites always appear in n-mer. A band of extended states 
will emerge in the binary alloy when the site energies are 
long-range correlated [23j. In addition, other theoretical 
models have been suggested to produce conducting states 
in low-dimensional disordered systems [24j-[33| , and some 
of them have been corroborated in GaAs-AlGaAs super- 
lattices [3] and Bose-Einstein condensate (BEC) [T3 |. 
Nevertheless, we notice that all electronic states become 
localized in these correlated disordered systems when the 
disorder degree is extremely large. 

In this Letter, we explore the electronic and localiza- 
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FIG. 1: (Color online) Schematic views of (a) random bi- 
nary alloy and of (b) RBL model. The former is composed 
of two different sites (atoms) A and B with constant nearest- 
neighbor hopping integral t, while the latter is constituted by 
assigning two types of layers La and Lb at random, where 
La (Lb) contains only A (B) sites. The coupling parameters 
within La and Lb layers are Xa and As, respectively, (c) 
Typical wavefunctions of a random binary chain (solid line) 
and of a corresponding binary layer with S = 3 (dashed line). 
The wavefunctions are obtained from direct diagonalization of 
the Hamiltonian with ea~0, £b=2, t=AA=As = l, and L=200 
[see Eq. ([1])] . All of the wavefunctions are bounded in a small 
region and both systems behave as insulators. However, a 
band of extended states will emerge in the RBLs by introduc- 
ing a relation between ea, A a and eb, As (see text). 



tion properties of a random binary layer (RBL) model, 
which can be constructed by coupling identical random 
binary chains along the transverse direction or equiva- 
lently by stochastically assigning two types of layers La 
and Lb in the lattice, as illustrated in Fig. [IJb). One 
expects that since a single random binary chain behaves 
as an insulator, all electronic states should be localized 
in the RBL model [see the dashed line in Fig. [He)]. How- 
ever, contrary to this physical intuition and the scaling 
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theory of localization, we show that a continuous band of 
extended states could emerge in the RBLs and the disor- 
dered systems behave as metals by introducing a relation 
for the model parameters, giving rise to an MIT in pa- 
rameter space. This is due to the presence of a perfectly 
ordered subband. Besides, the mobility edges can be de- 
termined analytically and the extended states are inde- 
pendent of the strengths of the diagonal and off-diagonal 
disorders. Our results still hold in two-dimensional (2D) 
and 3D disordered systems with arbitrary layer struc- 
tures, and may be implemented in the BECs. 

The Hamiltonian of the RBL model can be written in 
the tight-binding form: 
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where (cy) is the creation (annihilation) operator of 
an electron at site (i,j), with subscripts i E [1,5] labeling 
a chain and j 6 [1, L] denoting a layer. 5 is the number of 
chains and can also describe the cross section geometry, 
and L is the length. The site energy Ej for the La (Lb) 
layer is taken as ea (sb), the interchain hopping integral 
Xj is set to Xa (Xb) within the La (Lb) layer, and t is the 
intrachain hopping integral. (■ • • ) represents the nearest- 
neighbor sites. In what follows, the most disordered case 
is considered with probability A for each layer, ea = 
is chosen as the energy reference point, and all energy 
parameters are in units of t. Consequently, the diagonal 
disorder degree is W = £b- 

It is convenient to study the localization properties of 



the RBLs by using the transfer- matrix method [35 . 
i.e., the Lyapunov exponents j^s can be obtained from 
the QR decomposition together with the Gram-Schmidt 
reorthonormalization. We calculate the first 5 positive 
7i's, which refer to the 5 electron propagating modes 
along the longitudinal direction. The localization length 
£, measuring the extent of the wavefunction, is then de- 
fined as the reciprocal of the smallest positive tj. Our re- 
sults are averaged over an ensemble of disorder configura- 
tions to reach convergence. Finally, the zero-temperature 
conductance is G — (2e 2 / h) J2i = i l/cosh 2 (7iL) [3l |37|. 

In the site representation, the Schrodinger equation 
H\9) = E\9) can be expressed as 



(EI-H J )* J =t(* J _ 1 +* J+1 ). 



(2) 



Here, E is the Fermi energy, I is the 5 x 5 identity ma- 
trix, Wj is the sub-Hamiltonian matrix of the jth layer, 
and M/j = (if>ij, if>2ji • • • ,ipSj) T with rp^ the amplitude 
of the wavefunction at site (i,j) and T the transpose. 
It is clear that there are two different Hermitian matri- 
ces Hj's, namely Ha and Hb- Both Ha and Hb can 
be diagonalized by a single unitary matrix U through 
Pj = XJ^HjV, and the diagonal elements of Vj are the 
eigenvalues of Hj. Thus, Eq. ([2]) can be transformed into 



with $j = Ut*^ = (<f> 



Then the RBL 



model is decoupled into following Schrodinger equations: 



(E - Vkj)<t>ki = t(<f>kj-l + 4>kj+l), 



(4) 



where Vkj is the fcth eigenvalue of Hj. When the cross 
section is a line, the eigenvalue i/jy is 



= £j + 2Xj cos 



kn 
5 + 1' 



(5) 



when the cross section is a rectangle with S=S y x 5 2 , the 
Vkj with integer index k=(k y — 1)5 Z + k z is 



Vkj 
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Here, £j=0 and Xj=Xa when Hj=HA, £j—^B and Xj=Xs 
when Hj=H B , k E [1,5], k y e [1,5„], and k z E [1,5J. 
It clearly appears that the RBLs can be decoupled into 5 
chains, each of which is a random binary chain with site 
energies v^a and v^b determined by Eqs. (0 or ([6]); and 
the electronic and localization properties of the RBLs 
could be strongly changed by varying only, and the 
energy spectrum will be shifted. One notices that when 
vua — v kB , the fcth decoupled chain is a perfectly ordered 
one with constant site energy v^a and nearest-neighbor 
hopping integral t, and the electronic states are Bloch-like 
states in the energy interval [v^a ~ 2£j v^a + 2i], with the 
mobility edges being at E — v^a ± 2t. This result will be 
verified numerically for several RBLs with different cross 
sections (see below). 

We first consider the simplest case of 5=2, namely the 
two-leg ladder model [33|], which has been used to simu- 
late the charge transport through double-stranded DNA. 
In this respect, the site energies of two decoupled chains 
are vi A =Xa, v 1b =£b + X b and v 2 a=-Xa, v 2 b=£b-X b , 
respectively. Accordingly, the energy spectrum of the dis- 
ordered two-leg ladder, composed of the energy bands of 
individual decoupled chains, can be written as 



[via ~ 2£, via + 2t] U [v 1B - 2t, v 1B + 2t], 
[V2A - 2t, v 2A + 2t] U [v 2B - 2t, v 2B + 2i]. 
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Figures [21(a) and Efb) display respectively the density of 
states (DOS) and the localization length £ with £^=0.5, 
A,4=— 2.5, Xb=— 3 (solid lines), so that via — V\b- One 
can see from Fig. dja) that the energy spectrum con- 
sists of two subbands divided by an energy gap. The 
left subband resembles a periodic chain with extremely 
sharp van Hove singularities and smooth DOS profile, al- 
though both diagonal and off-diagonal disorders exist in 
the system; while the right subband refers to a disordered 
chain. This accords with Eq. {7J, indicating that each 
subband corresponds to an individual decoupled chain. 
On the other hand, it is clear that £ is comparable to 
L for the left subband [Fig. [21(b)] and the conductance 
G is exactly 2e 2 /h (see inset), implying a perfect prop- 
agating mode in this energy region. Actually, all states 
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FIG. 2: (Color online) (a) DOS and (b) localization length for 
disordered two-leg ladders with £s=0.5, Aa=— 2.5, Xb~— 3 
(solid lines), £b=0.5, Aa=Ab=— 2.5 (dashed lines) , and eb=0, 
Aa=— 2.5, A_b=— 3 (dotted lines). The inset displays the cor- 
responding conductance [38|. The DOS is calculated by diag- 
onalizing the Hamiltonian for a single disorder configuration 
and is almost the same when other disorder configurations are 
employed, since the system size is very large. In the gray en- 
ergy region, £ is averaged in a small energy window of width 
0.1 around E to avoid numerical fluctuations. 



in the left subband are extended in the thermodynamic 
limit, as verified by the scaling behaviors of the normal- 
ized localization length -| in Fig. (3][a) . We find that the 
scaling function /3=0 for all energies in the left subband 
Q and the disordered two-leg ladder exhibits metallic 
behavior. The states are localized elsewhere with G=0 
and the mobility edges are at E = Xa ± 2i. In fact, 
these results always hold when via = v\b or V 2A = V2B 
is satisfied, regardless of the diagonal and off-diagonal 
disorder degrees and the other parameters. Besides, the 
position of the energy band of extended states, relative to 
that of localized states, can be modulated by altering the 
model parameters [Eq. 0] and a variety of phenomena 
will occur, such as the shift of the mobility edges. These 
statements will be further discussed in the following. 

As a comparison, Fig.[5]shows DOS and £ by changing 
only one parameter of the two- leg ladder, viz., replacing 
Xb with —2.5 (dashed lines) and Eb with (dotted lines), 
respectively. In contrast, all states become localized for 
both two-leg ladders with G = (see inset) and con- 
stant localization length at fixed energies [from Fig. EJb) 
we obtain /3=— 1], due to the Anderson localization ef- 
fects, although they possess either the diagonal disorder 
(dashed lines) or the off-diagonal disorder (dotted lines) 
and are more ordered than the former case. 

Figures[3](c) and[3Jd) plot -| vs W (eb) for the two-leg 
ladders with Xa~— 2.5, Xb=— 2.5 — W and Aa=A_b=— 2.5, 
respectively. Contrary to the previous works that all 
states become localized in the disordered systems when 




FIG. 3: (Color online) Scaling behaviors of the normalized 
localization length at several energies for the corresponding 
two- leg ladders in Fig. [2] The states are extended in the ther- 
modynamic limit if * is invariant or increased by increasing L 
and are localized otherwise. £ vs L for (a) es=0.5, Xa~— 2.5, 
Xb~— 3 and for (b) es=0.5, Xa=Xb~~ 2.5. The dependence 
of j- on L can be fitted exactly by a simple function £ oc L 13 , 
as shown by the solid, dashed, and dotted lines, -g vs W (es) 
for (c) Xa=~2.5, A s =-2.5 - W and for (d) \ a =Xb=-2.5. 
The inset presents the enlarged view of j^-W in the region of 
W < W c with W C =E + 2t — A_g, and the corresponding fitting 
curves j oc . The parameter a in the vertical scale of (b) 
and (d) is employed to separate the lines, a— 2, 1, and 0.5 for 
£"=—2.5, —1.7, and —0.9, respectively. 



the disorder degree is very large 18|- 33 1 , we can see from 
Fig. [3jc) that 4 is independent of W for the former lad- 
der and all states are always extended in the gray energy 
region [Fig. |2fb)], because via = v\b=— 2.5 and the left 
subband denotes a completely ordered chain and does not 
move with W. While for the latter ladder, a crossover 
W c , dividing strong and weak dependence of |j on W, ex- 
ists in all the curves [Fig. EUd)], similar to that observed 
in the experiments [a Il3l . fl7j ] . ^ is strongly declined by 
increasing W up to W C =E + 2t — Xb, due to the gradu- 
ally enhanced scatterings from the potential barriers at B 
sites. The behavior of j vs W can be fitted by a power 
law relation for W < W c (see inset). 4 is slightly de- 
creased by further increasing W, since the energy E is 
separated from the band \v\b — 2i, v\b + 2t] for W > W c 
[see Eq. (J7}]. In this situation, the state at E is formed 
by A sites only and the electron is usually confined inside 
a single A site, leading to the localization length £ « 1. 

We then consider the RBLs of various cross sections. 
Figures Eff a )-S{d) show J vs E for several values of Eb , 
As, L by fixing S—5 and Xa=5. The eb and Xb are cho- 
sen so that VkA — VkB is satisfied for specific k and a per- 
fectly ordered chain always appears. A band of extended 
states as well as the mobility edges can be observed in 
all the situations and are shifted toward lower energies 
by increasing W; and the wavefunction in the metallic 
regime will spread over the entire system (see the dashed 
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FIG. 4: (Color online) Energy-dependent normalized localiza- 
tion length for the RBLs of linear and square cross sections 
with several values of L. 4- vs E for (a) eb=V3, Ab=4, for (b) 
£b=3, As=2, for (c) £b=— 4, As=l, and for (d) £b=— 5\/3, 
A s =0 with 5=5 and A A =5. ^ va E for (e) 5=2x2, e B =2, 
A A =3, A s =2 and for (f) 5=3x3, e s =2, A A =0, X B =V2. 
The inset presents the wavefunctions of two eigenmodes at 
£=-6.9901 (solid line) and £=-6.8110 (dashed line) in the 
metallic regime with L=200. The arrangement of A and B 
sites is the same as that in Fig. QJc) . 



line of the inset). In this regard, the RBLs may be used 
as switching devices and their electronic structure could 
be controlled by changing Eb and A b while keeping other 
parameters fixed. However, the localized state can sur- 
vive in the metallic regime (see the solid line of the inset), 
due to the superposition of extended and localized sub- 
bands, as can be inferred from Eq. ([5]). Similar results 
can be obtained for the RBLs of square cross sections, 
even when Aa=0 [Figs. Hte)-|3If)]. Therefore, we con- 
clude that a band of extended states will appear in the 
RBLs and the systems behave as metals by implement- 
ing a relation between ea, an d £_b, As, leading to an 
MIT in parameter space. These are irrespective of the 
diagonal and off-diagonal disorder strengths. 

Finally, we discuss the generalization of our results and 
the possible relevance to realistic systems, (i) In contrast 
to the random-dimer model that the extended states van- 
ish when it contains more than two different sites Tj^- 22 | , 
our results can still hold when the site energies are ran- 
domly selected from a certain probability distribution. 
For example, we consider Ej uniformly distributed within 
[— -y-, W-]. By properly choosing the interchain hopping 
integral that Vkj—E c always holds for specific k [Eq. ([3])] 



or k y , k z [Eq. ©] with E c a constant, a band of extended 
states will emerge in the energy region [E c — 2t, E c + 2t]. 
(ii) Our results remain valid for other 2D and 3D dis- 
ordered systems when v^j=E c is satisfied. For an SxL 
disordered system, the number of the condition Vkj=E c 
is 25 '~ 1 +(~ :L ) i except for the case cos -^-=0 which leads 
to an ordered system. Similarly for an SxSxL disor- 



dered system, the number is 



2S 2 -l+(-l) 



(iii) These 



results also hold for the disordered systems with various 
layer structures, such as rectangle and hexagon, as long 
as the ^fej's for fixed k (or k y , k z ) are constant or fol- 
low a periodic distribution by increasing j . (iv) Besides 
the double-stranded DNA, we believe our results to be 
relevant for the BECs, which offer an ideal platform for 
studying Anderson localization (l3-17|; and the model 
parameters can be precisely controlled |39{ . 

In summary, we suggest a scheme to produce metal- 
lic states in a series of disordered systems. Our results 
indicate that a band of extended states could emerge in 
the thermodynamic limit and the disordered systems be- 
have as metals by properly adjusting the model parame- 
ters, generating a metal-insulator transition in parameter 
space. The results are independent of the diagonal and 
off-diagonal disorder strengths, and still hold in two- and 
three-dimensional disordered systems. 
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